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Abstract 

We compute the analytic expression of the probability distributions Faex,+ 
and Faex- of the normalized positive and negative AEX (Netherlands) in- 
dex daily returns r(i). Furthermore, we define the a re-scaled AEX daily 
index positive returns r(t) a and negative returns (— r(t)) a that we call, after 
normalization, the a positive fluctuations and a negative fluctuations. We 
use the Kolmogorov-Smirnov statistical test, as a method, to find the values 
of a that optimize the data collapse of the histogram of the a fluctuations 
with the Bramwell-Holdsworth-Pinton (BHP) probability density function. 
The optimal parameters that we found are a + = 0.46 and a~ = 0.43. Since 
the BHP probability density function appears in several other dissimilar phe- 
nomena, our result reveals an universal feature of the stock exchange markets. 
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1. Introduction 



The modeling of the time series of stock prices is a main issue in economics 
and finance and it is of a vital importance in the management of large portfo- 
lios of stocks [13, Hal, H3 ■ 

Here, we analyze the AEX Index. The AEX index, 



derived from Amsterdam Exchange index, is a stock market index composed 
of Dutch companies that trade on Euronext Amsterdam, formerly known as 
the Amsterdam Stock Exchange. Started in 1983, the index is composed of a 
maximum of 25 of the most actively traded securities on the exchange. It is 
one of the main national indices of the pan-European stock exchange group 
Euronext alongside Brussels' BEL20, Paris's CAC 40 and Lisbon's PSI-20. 
The time series to investigate in our analysis is the AEX index from October 
of 1992 to October of 2009. Let Y(t) be the AEX index adjusted close value 
at day t. We define the AEX index daily return on day t by 

[) ~ Y(t-1) ' 

We define the a re-scaled AEX daily index positive returns r(t) a , for r(t) > 0, 
that we call, after normalization, the a positive fluctuations. We define the 
a re-scaled AEX daily index negative returns (—r(t)) a , for r(t) < 0, that we 
call, after normalization, the a negative fluctuations. We analyze, separately, 
the a positive and a negative daily fluctuations that can have different statis- 
tical and economic natures due, for instance, to the leverage effects (see, for 
example, 0,3 0,0)- Our aim is to find the values of a that optimize the 



data collapse of the histogram of the a positive and a negative fluctuations to 
the universal, non-parametric, Bramwell-Holdsofworth-Pinton (BHP) proba- 
bility density function. To do it, we apply the Kolmogorov-Smirnov statistic 
test to the null hypothesis claiming that the probability distribution of the a 
fluctuations is equal to the (BHP) distribution. We observe that the P values 
of the Kolmogorov-Smirnov test vary continuously with a. The highest P 
values P + = 0.59... and P = 0.31... of the Kolmogorov-Smirnov test are 
attained for the values a + = 0.46... and a~ = 0.43..., respectively, for the 
positive and negative fluctuations. Hence, the null hypothesis is not rejected 
for values of a in small neighborhoods of a + = 0.46... and oT = 0.43.... Then, 
we show the data collapse of the histograms of the a + positive fluctuations 
and a~ negative fluctuations to the BHP pdf. Using this data collapse, we do 
a change of variable that allow us to compute the analytic expressions of the 
probability density functions /aex + and Jaex - of the normalized positive 
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and negative AEX index daily returns 

fAEX,+ (x) = 4.39x- 054 W(20.28a; a46 -2.12) 
fAEX,-(x) = 3.52x- a57 W(17.53a; a43 -2.14) 

in terms of the BHP pdf Jbhp- We exhibit the data collapse of the histogram 
of the positive and negative returns to our proposed theoretical pdfs Jaex,+ 
and /aex-- Similar results are observed for some other stock indexes, prices 



of stocks, exchange rates and commodity prices (see [13j, ll4|). Since the BHP 
probability density function appears in several other dissimilar phenomena 
(see, for instance, 0, 0, S, EEI, 11, 16, 17, [2]]]), our result reveals an universal 



feature of the stock exchange markets. 

2. Positive AEX index daily returns 

Let T + be the set of all days t with positive returns, i.e. 

T + = {t: r(t) > 0}. 

Let n + = 2285 be the cardinal of the set T + . The a re-scaled AEX daily 
index positive returns are the returns r(t) a with t G T + . Since the total 
number of observed days is n = 4318, we obtain that n + /n = 0.53. The 
mean /i+ = 0.105... of the a re-scaled AEX daily index positive returns is 
given by 

ri^E r{t)a (1) 

ter+ 

The standard deviation cr+ = 0.049... of the a re-scaled AEX daily index 
positive returns is given by 



V E r ^ 2a - ^ 2 ( 2 ) 



n 

teT+ 

We define the a positive fluctuations by 

r+(t) = r(t) ";^ (3) 

for every t G T + . Hence, the a positive fluctuations are the normalized a re- 
scaled AEX daily index positive returns. Let = —1.97... be the smallest 
a positive fluctuation, i.e. 

L+ = min{r+(t)}. 

t€T+ 
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Let = 5.08... be the largest a positive fluctuation, i.e. 

i?+ = max{r+(t)}. 

teT+ 

We denote by F at+ the probability distribution of the a positive fluctuations. 
Let the truncated BHP probability distribution F B HP,a,+ be given by 



BHP,a,+ 



X 



Fbhp(x) 



Fbhp{FLu) — Fbhp(L+ 



where Fbhp is the BHP probability distribution. We apply the Kolmogorov- 
Smirnov statistic test to the null hypothesis claiming that the probability 
distributions F Qi+ and FBHP,a,+ are equal. The Kolmogorov-Smirnov P value 
P Qj+ is plotted in Figured! Hence, we observe that a + = 0.46... is the point 
where the P value P a+ = 0.59... attains its maximum. 



Figure 1: The Kolmogorov-Smirnov P value P a ,+ for values of a in the range [0.30, 0.60]. 



It is well-known that the Kolmogorov-Smirnov P value P a ,+ decreases with 
the distance \\F at+ — Fbhp,o,+ \\ between F a ^ + and FsHP,a,+- In Figure |2l we 
plot D a + t+ (x) = \F a + )+ (x) — F B HP,a+,+ (%) I and we observe that D a + + (x) 
attains its highest values for the a + positive fluctuations above or close to 
the mean of the probability distribution. 
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Figure 2: The map D 0A6i+ (x) = \F 0A6 . + (x) - F B hp,qa6,+( x )\ 



In Figures |3] and HI we show the data collapse of the histogram f a +,+ of the 
a + positive fluctuations to the truncated BHP pdf fBHPa+ +• 




Figure 3: The histogram of the a + positive fluctuations with the truncated BHP pdf 
fBHP,QA6,+ on top, in the semi-log scale. 



Assume that the probability distribution of the a + positive fluctuations 
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Figure 4: The histogram of the a + positive fluctuations with the truncated BHP pdf 

fBHP,0A6,+ On top. 



r+ + (t) is given by F BHPja+j+ (see [nj). The pdf f A EX,+ of the AEX daily 
index positive returns r(t) is given by 

_ a+x<* + ~if BHP ((a- + - n+ + ) /a+) 

fAExAx) = PwftHwl^ ' 

Hence, taking ot + = 0.46..., we get 

Iakx ■[■'■) = 4.39...X- - 54 - W(20.28...x - 46 - - 2.12...). 

In Figures [5] and O we show the data collapse of the histogram / 1+ of the 
positive returns to our proposed theoretical pdf Jaex,+- 

3. Negative AEX index daily returns 

Let T~ be the set of all days t with negative returns, i.e. 

T~ = {t : r(t) < 0}. 

Let rr = 2021 be the cardinal of the set T~. Since the total number of 
observed days is n = 4318, we obtain that n~/n = 0.47. The a re-scaled 
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Figure 5: The histogram of the fluctuations of the positive returns with the pdf Jaex,+ 
on top, in the semi-log scale. 




Figure 6: The histogram of the fluctuations of the positive returns with the pdf Jaex,+ 
on top. 
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AEX daily index negative returns are the returns (— r(t)) a with t G T~. We 
note that — r(t) is positive. The mean fi~ = 0.122... of the a re-scaled AEX 
daily index negative returns is given by 

A£ = ^£(- r (*)r (4) 

teT- 

The standard deviation a~ = 0.057... of the a re-scaled AEX daily index 
negative returns is given by 



-1 (- r (t)) 2a ~ (^) 2 (5) 



n 

teT 



We define the a negative fluctuations by 

r _ (t) = Hte£ ( 6) 

for every t G T _ . Hence, the a negative fluctuations are the normalized a re- 
scaled AEX daily index negative returns. Let L~ = —1.97... be the smallest 
a negative fluctuation, i.e. 

L~ = min{r-(t)}. 

teT- 

Let R~ = 4.13... be the largest a negative fluctuation, i.e. 

R~ = max{r-(t)}. 

We denote by F a _ the probability distribution of the a negative fluctuations. 
Let the truncated BHP probability distribution Fbhp^,- be given by 

F (r) F BHP {x) 



FBHp{R a ) — FBHp{L a 



where Fbhp is the BHP probability distribution. We apply the Kolmogorov- 
Smirnov statistic test to the null hypothesis claiming that the probability 
distributions F a _ and F B HP,a- are equal. The Kolmogorov-Smirnov P value 
P a ,- is plotted in Figured Hence, we observe that a~ = 0.43... is the point 
where the P value P a - = 0.31... attains its maximum. The Kolmogorov- 
Smirnov P value P a _ decreases with the distance \\F a _ — F B HPa -II between 
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Figure 7: The Kolmogorov-Smirnov P value P a - for values of a in the range [0.30, 0.60]. 
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Figure 8: The map D 0A3 _(x) = |F . 4 3,-(^) - F B hp,oa3,-(x)\. 
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F a _ and F BHPa _. InFigure[SJ we plot D a -_(x) = \F a -_(x) — F B Hp, a -,-( x )\ 
and we observe that D a - _(x) attains its highest values for the a~ negative 
fluctuations below the mean of the probability distribution. 

In Figures 191 and [TO] we show the data collapse of the histogram /„- _ of the 
a" negative fluctuations to the truncated BHP pdf f B HP,a-,-- 




Figure 9: The histogram of the a negative fluctuations with the truncated BHP pdf 
fBHP,o.43 - on top, in the semi- log scale. 



Assume that the probability distribution of the a negative fluctuations 
r~_(t) is given by F BHP)0l -_, (see The pdf Jaex- of the AEX daily 



index (symmetric) negative returns —r(t), with T G T , is given by 



f. 



AEX, 



(X) 



a~x a l f B HP | 1 • 



a--{F BHP (R-_)-F BHP (L-_))- 
Hence, taking a~ = 0.43..., we get 

f A EX-{x, 



3.52...x'°- 57 -f BH p(l7.53...x 0A3 - - 2.14... 



In Figures [TT] and [121 we show the data collapse of the histogram fi of the 
negative returns to our proposed theoretical pdf Jaex-- 
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Figure 10: The histogram of the a negative fluctuations with the truncated BHP pdf 

fBHP,0.43- On top. 





Figure 12: The histogram of the negative returns with the pdf Jaex,- on top, in the 
semi-log scale. 



4. Conclusions 

We used the Kolmogorov-Smirnov statistical test to compare the his- 
togram of the a positive fluctuations and a negative fluctuations with the 
universal, non-parametric, Bramwell-Holdsworth-Pinton (BHP) probability 
distribution. We found that the parameters a + = 0.46... and a~ = 0.43... for 
the positive and negative fluctuations, respectively, optimize the P value of 
the Kolmogorov-Smirnov test. We obtained that the respective P values of 
the Kolmogorov-Smirnov statistical test are P + = 0.59... and P~ = 0.31.... 
Hence, the null hypothesis was not rejected. The fact that a + is different 
from a~ can be do to leverage effects. We presented the data collapse of 
the corresponding fluctuations histograms to the BHP pdf. Furthermore, we 
computed the analytic expression of the probability distributions Faex,+ an d 
Faex- of the normalized AEX index daily positive and negative returns in 
terms of the BHP pdf. We showed the data collapse of the histogram of the 
positive and negative returns to our proposed theoretical pdfs Jaex,+ and 
Jaex-- The results obtained in daily returns also apply to other periodici- 
ties, such as weekly and monthly returns as well as intraday values. 

In 13, [9J, it is found the data collapses of the histograms of some other 



stock indexes, prices of stocks, exchange rates, commodity prices and energy 
sources 



121 to the BHP pdf. 



Bramwell, Holdsworth and Pinton [3] found the probability distribution 
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of the fluctuations of the total magnetization, in the strong coupling (low 
temperature) regime, for a two-dimensional spin model (2dXY) using the 
spin wave approximation. From a statistical physics point of view, one can 
think that the stock prices form a non-equilibrium system 0, 18, 19, 23 



Hence, the results presented here lead to a construction of a new qualitative 
and quantitative econophysics model for the stock market based in the two- 
dimensional spin model (2dXY) at criticality (see fl~3|). 
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Appendix A. Definition of the Bramwell-Holdsworth-Pinton prob- 
ability distribution 

The universal nonparametric BHP pdf was discovered by Bramwell, Holdsworth 
and Pinton j^]. The BHP probability density function (pdf) is given by 



fBHp(v) 



dx 
2^ 



\ 



1 v-^ 1 lx V 

2^ e 



2N2 tl A * 



E 



N-l ix_ J_ 
k=l 2N At 



-farctanl^f- 



v^JV-l 1 



iln( i+-S^- 



(Al) 



where the {Afc}^ =1 are the eigenvalues, as determined in [5j, of the adjacency 
matrix. It follows, from the formula of the BHP pdf, that the asymptotic 
values for large deviations, below and above the mean, are exponential and 
double exponential, respectively (in this article, we use the approximation 
of the BHP pdf obtained by taking L = 10 and N = L 2 in equation ( lA.lj) ). 
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As we can see, the BHP distribution does not have any parameter (except 
the mean that is normalize to and the standard deviation that is normal- 
ized to 1) and it is universal, in the sense that appears in several physical 
phenomena. For instance, the universal nonparametric BHP distribution is 
a good model to explain the fluctuations of order parameters in theoretical 
examples such as, models of self-organized criticality, equilibrium critical be- 
havior, percolation phenomena (seejij), the Sneppen model (see and Q), 



and auto- ignition fire models (see [24j). The universal nonparametric BHP 



distribution is, also, an explanatory model for fluctuations of several phe- 
nomenon such as, width power in steady state systems (see fluctuations 
in river heights and flow (see 0, 0, [l5|, [ll|, E]), for the plasma density fluc- 
tuations and electrostatic turbulent fluxes measured at the scrape-off layer 
of the Alcator C-mod Tokamaks (see [25J) and for Wolf's sunspot numbers 



fluctuations (see |l7l|). 
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